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OPTIMIZATION OF THE MANAGEMENT OF 
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‘igure 1-1 : Vallée de la Vesdre du barrage d'Eupen à la confluence avec l'Ourthe et la Meuse (Dewals, et 
al., 2011) 


INTRODUCTION 


The effects of climate change are as manifold and various than the 
variation of the oceans level, the migration of some kind of species,... 
The goal of this work is to study the effects of climate change on 


dams management. 
This study is applied on two structures located in the Vesdre basin, in 


Belgium. From an integrated model including a climate model, a 
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hydrological model of the Vesdre watershed, a hydraulic calculation 
model of this watercourse and a model of dams management, the 
Vesdre discharges and the lakes levels can be calculated. Thanks to 
climate data feeding the model, a study of the influence of some 
parameters of the dam management rule is done, after having 
defined some original performance indicators. These indicators 
represent the damages due to floods, the hydropower production, 
the support of low flows and the production of drinking water 
throughout the year. Afterwards, from two climate scenarios, a dry 
and a wet, the changes induced in the Vesdre valley and in the tanks 
water levels are determined at the time horizons 2020-2050 and 
2070-2100. Then, the sensitivity analysis, performed in the period 
1974-2004, is used to research adaptations in the management rule 
to mitigate the effects of climate change 
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Figure 1-2 : Représentation schématique du réseau hydrologique de la Vesdre, de sa source jusqu'à sa 
confluence avec la Hóegne 


1-The multiple functions of a tank 


During the construction of a water retention structure, the objectives 
which justified this construction can be many and varied. These 
different objectives are conventionally grouped into two main 
functions: reducing the impacts of downstream flooding and 
mitigating the consequences of a dry season on water volumes 
(Dewals, et al., 2011). 

First, during flooding, the water runoff on the catchment area 
relating to the dam can be stored by the latter, provided that it has 
sufficient retention capacity called a sludge reserve. Thus, by 


temporarily preventing this volume of water from reaching the river 
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downstream, the hydrograph relating to the watershed whose 
waters are retained by the dam can be subtracted from the total 
flood hydrograph, in the absence of the ‘structure at any point 
downstream of the reservoir. This is flood clipping (Figure 2-1). Once 
the flood has passed, it is necessary to restore the sludge reserve by 
releasing water downstream, at a time when the danger of flooding 
has been eliminated. Obviously, this clipping of the flood will be all 
the more efficient as the watershed relative to the reservoir will 


represent a significant part of the entire watershed of the 











watercourse. 
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Figure 2-1: Illustration de la fonction d’écrétage de crue des barrages de la Vesdre sur l'hydrogramme de 
crue en aval des ouvrages de rétention 
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Indeed, the objective of maximum filling of the 
reserve, in anticipation of the dry season, is 
opposed to the objective of creating a sludge 


reserve Ime C ess E IO eio oae sau Seed 





compromise which consists in ensuring a sufficient 
reserve for the storage of the flood that must be 
capped, by a possible preventive operation of 
emptying the backwater, without jeopardizing the 
objective of supplying water in dry periods. 
Likewise, the low flow support function is opposed to 
the drinking warer production funetron Following che 
water loss it causes. 

Figure 2-2 clearly shows the main variables involved 
in optimizing the management of a reservoir. 
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Figure 2-2 : Variables et valeurs seuils influencant les opérations de gestion d'un ouvrage de retenue d'eau 
(Bieri, et al., 2013) 
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2-Methods of optimizing reservoir 


management 


Theifrurturc hydrological conbtrrbusronms ob che 





watershed relating to the studied reservoirs are 


generally used to define an appropriate management of 


these reservoirs. By their nature, these inputs are 


random and can only be predicted with some precision 


in the very near future. 


Thus, optimal management of 


reservoirs is complex to define since it is based on 


a Ppaertetal agnorance or the Short-term input 


parameters and on a significant ignorance of these 


long-term parameters. 
The study of the optimization of 
the production of hydroelectric 
power was approached in 1923 
by Varlet, graphically by the 
Tense Wire method. In 1946, 
Masse provided proof of the 
effectiveness of the method 
proposed by Varlet through an 
economic and mathematical 
analysis of the marginal benefit 


derived from the production of 


hydroelectric power ona simple 
case. This Massé model will be 
generalized to several reservoirs 
and to several objectives 
thereafter, in a limited future 
(Morel-Seytoux, 1998). A review 
of the state of the art 
concerning the optimal 
management of reservoirs was 
carried out by J. W. Labadie 
(2004). The problem arises from 
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an objective function 
incorporating criteria that may 
be economic, environmental, 
ecological or social. Besides this 
objective function, various 
constraints are defined relating 
to the extreme levels of the 
lakes, the low flow rates 
downstream of the reservoirs, 
etc. When possible, a definition 
of the objective function in 
terms of costs simplifies the 
resolution. of the problem by 
finding the minimum of this 
function (Watkins, et al., 2006). 
Various cost functions were 
developed for this purpose 
(Jordan, 2007). There is a large 
number of methods allowing 


the resolution of the problem 


posed: implicit stochastic 
optimization, linear and 
nonlinear programming models, 
stochastic dynamics,... These 
various models are also 
presented in the article by JW 
Labadie (2004), but will not be 
not used for this job. The 
concern is that management 
objectives cannot always be 
quantified in monetary units. 
Consider the environmental 
consequences of low flow rates 
or the social consequences of a 
shortage of drinking water 
supply. In this case, the 
definition of an objective 
function to be minimized can be 
replaced by a multicriteria 


analysis 


3-Local sensitivity analysis 
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The presentation of local sensitivity analysis, as 
used in this work, is mainly based on the thesis of 
S. Wildemeersch (2012). It is used in the relaxation 
phase of the parameters in order to determine their 
respective influences on the performance indicators. 
Local sensitivity analysis involves studying a system 
around a limited number of parameters studied 
independently. Theoretically, this method is only 
valid an linear models, but practice has shown i: toe 
be acceptable when the model is not "extremely non- 
linear" (Hill, et als, 2007). The advantage of this 
method is that it considerably reduces the number of 
Simulations to be carried out. However, the order of 
magnitude of a simulation being several days, this 
Simplurucakvom vseressential am tie contexu OF hrs 
work 

After having chosen a certain number of parameters bi 
to vary and the performance indicators yj 
representative of the tank objectives, the 
sensitivity matrix J can be created by a number of 
simulations Corresponding Eto ine number or 


parameters: 
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aby|, ODrpar |, 
J= : 3 : 
OYnind| OYnind 


With npar the number of parameters and nind the number 
ONE IMCL ecu eme 

Once the matrix J is completed, the influences of the 
parameters on the indicators are analyzed by the dss 
for dimensionless scaled sensitivities and by the css for 
composite scaled sensitivities. The elements dssij 
assess the importance of the variations of indicator 
i as a function of parameter j, by dimensionless 
formatting. The elements cssij are a measure of the 
Sensitivity of all the indicators to a parameter j 
and therefore make it possible to choose the most 


influential parameters to keep for optimization. 





dy; 1 
dss, = | xlb|x | (2-2) 
y ab A | | yi 
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(x (dssj |) 


1 - nind 


(2-3) 
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4-Integrated model 
basin 


Optimizing the management of 


the Vesdre dams, with a view 





to reducing the damage caused 
by floods, guaranteeing a 
continuous supply of water 
throughout the year, 
producing a maximum quantity 
of hydroelectric energy and 
supporting the low flows, 
taking into account climate 
change, requires a number of 
numerical models listed below 
(OL Cie ells, 2000o 

e A climate model for 
determining precipitation and 


temperature; 


e A hydrological model for 





converting climatic data into 


inflow flows over the entire 





watershed of the watercourse; 


of the Vesdre 


e A reservoir management 


model 





(developed in the next 
chapter); 

e A hydraulic model of the 
watercourse with, as inputs, 
the inflow rates of the 
hydrological model; 

e Possibly, a flood risk 


calculation model allows a 


more precise quantification 





of the objective of limiting 


flood damage. 





All the climatic, 
hydrologruealy, saydraul vey, 
reservoir management and 
flood damage calculation 
models are grouped into an 
integrated model of the 
Vesdre basin, developed 
within the framework of the 


AMICE project (Detrembleur, 
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analyzes of climate change 
and the sensitivity of the 
parameters of the management 


law of the Eupen and Gileppe 


dams do not rcguirs the rull 








use of the integrated model 
of the Vesdre basin 
developed. Indeed, based on 
the results provided by the 
integrated model of the AMICE 
project, the use of a one- 
dimensional hydraulic model, 
fed by the climatic and 
hydrological models, is 


sufficient to update the 


3-1 Climate model 


The presentation of the 
climate model is based 
on a report from the 
AMICE project (Drogue, 


Cb al. ZOOR EAS 


results of the complete 











model. It goes without saying 





hate they gain Inm eSrmS Of 








computation time is 





substantial and justifies 
this approach. The different 


models constituting the 





integrated model of the 
Vesdre basin are first 
presented below. Then, the 


process of updating the data 





necessary for studies of the 
influences of the parameters 
of the management law and 
analyzes of the effects of 


climate change is presented. 


model uses the Delta 
change approach, 
considered robust and 
ELANSDarent, tO account 


for climate change (te 
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made et al 201i 
korci, et akl, 2007) 
ibis merkod  consises im 
projecting the 
precipitation and 
temperatures, measured 
in recent years, to the 
Ere opos ED) 
and 2090-2100 from 
seasonal trends 
influenced by the 


greenhouse gas emission 


GCM/RCM 


SRES 
Emission Sce- > outputs E 


narios (daily time step) 





Scenarios selected, 
among the set of 
Scenarios developed by 
che ppc X DEO ZOO) 
(Figure 3-1). The use of 
regional climate models 
(RCM) allows to refine 
the results obtained by 
global climate models 
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Figure 3-1 : Génération de scénarios climatiques à partir de la méthode du changement delta (Drogue, et 


All the scenarios 
developed by the IPCC 
being equiprobable at 
Present Ewo extreme 


scenarios, among the 


multitude of scenarios 
Simulated by tene TELE, 
one dry and one wet, 


were chosen within the 
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framework of the AMICE 
project. 
On the basis of these 
scenarios, Projections 
concerning precipitation 
and temperatures, in the 
Walloon part of the 
Meuse watershed, were 
made for the two time 


horrzOonsSes Etude den 
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sigress tc isearo NE BE 
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of the greenhouse gas 
emission scenario 
chosen. Then, 
precipitation is reduced 
zn summer in both 
Scenarios and increased 
in winter in the wet 


Scenario. 
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Figure 3-2 : Tendances saisonnières des précipitations et températures, selon le scénario utilisé, et pour 
les horizons temporels 2020-2050 (foncé) et 2070-2100 (clair) (Drogue, et al., 2010) 


To illustrate the principle of Delta change, Figure 


323 compares the flows at Chaudtontaine, during Ehe 
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major flood of September 1998, with the humid climate 
scenario for the 2070-2100 horizon, created by 
ALSrupLion Of current iclimare data. 
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igure 3-3 : Principe du changement Delta illustré sur les débits à Chaudfontaine en septembre 1998 


3-2 : Description and modeling of the 
current management procedure 


It is important to stress that the current management 
Odi damsecamepEdol aqace Mie nor Automatic, bur restes 
from decisions made by one or more people, however 


based on the procedure described in the handling 
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note. The simulation created therefore inevitably 





deviates from the actual management 





Figure 4-1 : Résumé dp la modélisation initiale de la loi de gestion des barrages de la Vesdre. 
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